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Abstract. The bulk viscosity of the neutron star mat- 
ter due to the direct Urea processes involving nucleons, 
electrons and muons is studied taking into account possi- 
ble superfluidity of nucleons in the neutron star cores. The 
cases of singlet-state pairing or triplet-state pairing (with- 
out and with nodes of the superfluid gap at the Fermi sur- 
face) of nucleons are considered. It is shown that the super- 
fluidity may strongly reduce the bulk viscosity. The prac- 
tical expressions for the superfluid reduction factors are 
obtained. For illustration, the bulk viscosity is calculated 
for two models of dense matter composed of neutrons, 
protons, electrons and muons. The presence of muons af- 
fects the bulk viscosity due to the direct Urea reactions 
involving electrons and produces additional comparable 
contribution due to the direct Urea reactions involving 
muons. The results can be useful for studying damping of 
vibrations of neutron stars with superfluid cores. 



1. Introduction 

The dissipative processes in neutron stars play an impor- 
tant role for some dynamical properties of these unique 
objects. Shear viscosity damps differential rotation of neu- 
tron stars, leading to their uniform rigid-body rotation. 
Quite generally, the viscosity of neutron star matter im- 
plies damping of pulsations of neutron stars. Such pulsa- 
tions could be excited during the process of neutron star 
formation. They could also be triggered by instabilities ap- 
pearing during neutron star evolution, or by some external 
perturbations. The corresponding damping timescales in- 
volve the shear and bulk viscosities of neutron star in- 
terior. Both viscosities depend on density, temperature 
and composition of dense matter. Calculations of damping 
timescales of pulsations for various models of neutron star 
interiors have been done by Cutler et al. ( 199C| ). Viscous 
damping of pulsations of newly born hot neutron stars 
turns out to be due to the bulk viscosity. 

Another role of the viscosity of neutron star matter 
is that it can damp gravitational radiation driven insta- 
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bilities in rotating neutron stars and, therefore, could be 
important for determination of the maximum rotation fre- 
quency of neutron stars. In the absence of viscosity all ro- 
tating neutron stars would be driven unstable by the emis- 
sion of gravitational waves. Viscous damping timescales 
enter explicitly the criteria for the appearance of these in- 
stabilities. Similarly, as in pulsating non-rotating neutron 
stars, viscous damping of gravitational radiation driven 
instabilities in rapidly rotating newly born neutron stars 
is dominated by bulk viscosity of neutron star interiors 
(e.g., Lindblom |l99|, Zdunik 1996, Lindblom et al. |1998D . 

In this paper, we focus on the viscosity of matter in 
the neutron star cores (which extend from the layers of 
density p ~ 1.5 x 10^'' g cm~^ to the stellar centers). It is 
well known that the cores consist of baryons (neutrons n, 
protons p and possibly hyperons) and leptons (electrons e 
and muons fi) . All constituents of matter are strongly de- 
generate fermions. The electrons and muons form almost 
ideal Fermi gases. The electrons are ultrarelativistic while 
the muons may be non-relativistic or relativistic depend- 
ing on density. The nucleons are, to a good approxima- 
tion, non-relativistic and constitute strongly interacting 
Fermi liquid. At the densities close to the normal nuclear 
density (baryon number density hq = 0.16 fm^"^ which 
corresponds to the mass density po = 2.8 x 10^^ g cm^'^), 
neutron star matter is composed of n, p, e, and fj,. At 
higher densities [p ^ (3-4) po] some models of dense mat- 
ter predict appearance of hyperons. At still higher den- 
sities, some calculations indicate possible presence of ex- 
otic phases (pion condensate, kaon condensate, deconfined 
quark matter). We will not consider the hyperonic or ex- 
otic phases but restrict ourselves to the study of the npefi 
matter. 

Our analysis is additionally complicated by possible 
superfluidity of nucleons in the neutron star cores. The 
superfluidity is thought to be produced by Cooper pair- 
ing of nucleons due to attractive parts of nucleon-nucleon 
interaction. The superfluidity of nucleons in the neutron 
star cores has been the subject of numerous papers (as 
reviewed, for instance, by Yakovlev et al. 1999). Various 



microscopic theories predict very different superfluid gaps 
(critical temperatures T^n and T^p) of neutrons and pro- 
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tons depending on specific model of strong interaction em- 
ployed and specific many-body theory used to account for 
medium effects. However, all these results have important 
common features. In particular, the proton pairing occurs 
mainly in the ^So-state since the pp interaction is attrac- 
tive in this state everywhere in the neutron star core due 
to not too high number density of protons. As for the neu- 
trons, their interaction in the ^Sq state turns from attrac- 
tion to repulsion at densities p ^ po but the interaction in 
the '^P2 state may be attractive and may lead to Cooper 
pairing. The critical temperatures T^n and Tcp in the neu- 
tron star cores predicted by different microscopic theories 
depend on density and scatter in a wide range from about 
10® to 10^" K. Under these conditions we will not rely on 
any specific microscopic theory of nuclcon superfluidity, 
but will treat T^n and Tcp as free parameters. 

The viscosity, we are interested in, is well known to 
consist of the shear viscosity and bulk viscosity. The stan- 
dard source of the shear viscosity of the neutron star mat- 
ter is scattering between its constituents. Classical calcula- 
tions of shear viscosity for the npe model of non-superfluid 
matter were done by Flowers & Itoh (1979). Their results 
were used in the studies of damping of neutron star pul- 
sations by Cutler et al. (199C). In the superfluid core of a 
rotating neutron star, there is an additional viscous mech- 
anism, called mutual friction, resulting from the scattering 
of electrons off the magnetic field trapped in the cores of 



superfluid neutron vortices (Lindblom & Mendell 1995). 

The bulk viscosity may partly be determined by parti- 
cle scattering. However, this component of bulk viscosity is 
usually m uch sm aller than the shear viscosity (e.g., Baym 
& Pethick 1991 ) . The main contribution into the bulk vis- 
cosity of sufficiently hot npep matter comes from the neu- 
trino processes of Urea type associated with electron and 
muon emission and capture by nucleons. We will focus on 
such bulk viscosity. Generally, the neutrino processes in 
question are divided into the direct Urea and the modified 
Urea processes. A direct Urea process is a sequence of two 
reactions (direct and inverse one) and can be written as 



p + I + I'l, p + I ^ n + vi, 



(1) 



where / is either electron or muon, and vi is an associated 
neutrino. The most important is the process (Lattimer et 
al. 1991) involving electrons (Z = e); it consist of the beta 
decay of neutron and subsequent beta capture. It should 
be emphasized that the both direct Urea processes are 
forbidden by momentum conservation of reacting parti- 
cles for the simplest model of dense matter as a gas of 
noninteracting Fermi particles (e.g., Shapiro & Teukolsky 
198^ ) at any density p in the neutron star cores. Never- 



theless they are allowed (Lattimer et al. 1991) for some 
realistic equations of state at densities higher than some 
threshold densities (of several po)- Thus, the direct Urea 
processes may be open in the inner cores of rather massive 
neutron stars. The threshold density for the muon process 
is always higher than for the electron one. 



If allowed, the direct Urea processes produce the most 
powerful neutrino emission from the neutron star cores 
(Lattimer et al. |199l| ). Corresponding neutrino emissivi- 
ties were calculated by Lattimer et al. ( |1991| ) and used 
in numerous simulations of the neutron star cooling as 
reviewed, for instance, by Yakovlev et al. ( [L999[ ). In the 
absence of nucleon superfluidity, the direct Urea processes 
lead to the fast cooling of neutron stars. If allowed, the 
direct Urea processes produce the main contribution into 
the bulk viscosity we are interested in. 

However, for many equations of state the direct Urea 
processes are forbidden by momentum conservation at any 
density in the neutron star cores. Moreover, they are pro- 
hibited at p 3 po for the majority of equations of state. 
In such cases, they do not operate in the low and medium- 
mass neutron stars and in the outer cores of all neutron 
star models constructed using these equations of state. If 
so, the bulk viscosity is determined by the the reactions 
of the modified Urea processes 

n + N^p + N + l + Di, p + N + l^n + N + vi, (2) 

where N is an additional nucleon required to conserve 
momentum of the reacting particles. For instance, in npe 
matter one has two modified Urea processes corresponding 
to N — n and N — p, respectively, which can be referred 
to as the neutron and proton branches of the modified 
Urea process (e.g., Friman & Maxwell 1979| , Yakovlev & 
Levenfish 1995). The rates of the modified Urea processes 



are typically 3-5 orders of magnitude lower than the rates 
of the direct Urea processes. The modified Urea processes 
either have no density threshold (as the neutron branch in 
npe matter) or have much lower density thresholds than 
the direct Urea processes. Thus they operate in the entire 
neutron star core. If the direct Urea processes are for- 
bidden and matter is non-superfluid, the modifled Urea 
processes produce the main neutrino emission from the 
neutron star cores leading to slow {standard) cooling of 
neutron stars. Their role in the neutron star cooling the- 
ory h as be en studied in many papers (see, e.g., Yakovlev 
et al. 199£ , for review). 

Thus, the problem of calculating the bulk viscosity 
due to neutrino processes is quite complicated: there are 
several neutrino processes involved influenced by possible 
nucleon superfluidity. So far the bulk viscosity has been 
studied only for non-superfluid npe matter. The viscosity 
due to the neutron branch of the modified Urea process 



was analyzed by Sawyer (1989) while the viscosity pro- 
duced by the nucleon direct Urea process was considered 



by Haensel & Schaeffer (1992). The effects of superffuidity 
have not been analyzed for the problem of bulk viscosity 
but studied thoroughly for the neutrino emissivity pro- 
duced in different reactions (e.g., Yakovlev et al. 1999, 
and references therein). 

The relative importance of the bulk viscosity produced 
by neutrino reactions with respect to the shear viscosity 
produced by collisions can be estimated by comparing the 
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results by Sawyer ( |l98g| ) and Haensel & Schaeffer ( |1992D 
with the values of the shear viscosity calculated by Flow- 
ers & Itoh ( p^79| ). The comparison shows that the neu- 
trino bulk viscosity dominates in the neutron star cores for 
temperatures T <; 10^ K if the direct Urea processes are 
switched on and for T ^ 10^ K if the direct Urea processes 
are forbidden. In superfluid matter, the bulk viscosity can 
be even more important. 

In this paper, we consider the bulk viscosity produced 
by the direct Urea processes in npe^ matter of the neu- 
tron star cores. In analogy with the effects of superfluidity 
on the neutrino emissivity, we will analyze the effects of 
superfluidity of nucleons on the bulk viscosity. 



2. Bulk viscosity in non-superfluid matter 
2.1. Bulk viscosity in npefi matter 

Consider the bulk viscosity produced by the direct Urea 
process (involving muons and electrons) in non-superfluid 
npefi matter. 

Due to very frequent collisions between particles, dense 
stellar matter should very quickly (instantaneously on 
macroscopic time scales) achieve a quasi-equilibrium state 
with certain temperature T and chemical potentials fii of 
different particle species i = n,p, e,/j,. Typically, all par- 
ticle species are strongly degenerate. We assume that the 
matter is transparent for neutrinos, which therefore do not 
contribute to the thermodynamical quantities. 

A quasi-equilibrium state described above does not 
mean full thermodynamic equilibrium. The latter assumes 
additionally the equilibrium with respect to the beta and 
muon decay and capture processes. We will call it the 
chemical equilibrium. Relaxation to the chemical equilib- 
rium depends drastically on a given equation of state and 
local density of matter p. It is realized either through di- 
rect Urea or modifled Urea processes (Sect. 1). Consid- 
eration of this subsection is valid for all Urea processes 
although the practical expressions (Sects. 2.3 - 2.4) will 
be obtained for the direct Urea processes. 

The Urea processes of both types, direct and modified, 
are rather slow. Although the chemical relaxation rate de- 
pends strongly on temperature, in any case it takes much 
more time (from tens of seconds to much longer time in- 
tervals) than the rapid relaxation to a quasi-equilibrium 
state described above. Therefore, a neutron star can be in 
a quasi-equilibrium, but not in the chemical equilibrium, 
for a long time. 

If the chemical equilibrium is achieved, then the chem- 
ical potentials satisfy the equalities /i„ — j-ip -\- fie and 
ti'n — fJ'p + fJ'n, which imply fie — Mm- Under these condi- 
tions, the rates of the direct and inverse reactions, F/ and 
F/ (Z = e or /x), of any Urea process are equal. 

Let us assume that the neutron star undergoes radial 
pulsations of frequency lu. Associated temporal variations 
of the local baryon number density will be taken in the 



form rifc — n^o + ribicosujt, where nbi is the pulsation 
amplitude and rtbo is the non-perturbed baryon number 
density (|rt6i| <C ?i6o)- We assume further that rtbo corre- 
sponds to the chemical equilibrium. This chemical equi- 
librium is violated slightly in pulsating matter. If the pul- 
sation frequency lu were much smaller than the chemical 
relaxation rates, the composition of matter would follow 
instantaneous values of Ub, realizing the chemical equilib- 
rium every moment of time. 

In reality, the typical frequencies of the fundamental 
mode of the radial pulsations uj ~ lO'^-lO'* s~^, are much 
higher than the chemical relaxation rates. As a result, 
the partial fractions Xi — Ui/rib of all the constituents 
of dense matter are almost unaffected by pulsations (i.e., 
almost constant). Owing to the slowness of the Urea reac- 
tions, these fractions lag behind their instantaneous equi- 
librium values, producing non-zero differences of instan- 
taneous fii: 



f^n f^p f^ei ^/j f-^n f^p M/i- 



(3) 



This causes an asymmetry of the direct and inverse di- 
rect Urea reactions, and, hence, slight deviations from the 
chemical equilibrium. The asymmetry, calculated in the 
linear approximation with respect to rji , is given by 



Ti-Ti = -A; rfi 



(4) 



where A/ are the coefficients specified in Sect. 2.4 for the 
direct Urea reactions. Microscopic calculation (Sect. 2.4) 
yields Ag = A^. In this paper, we will restrict ourselves to 
the case \rfi\ <^ T {I = e, fi). Our d efinition of A/ is the 



same as was used by Sawyer ( 1989 ) for the case of npe 
matter. Notice that A; defined in this way is negative. 

The non-equilibrium Urea reactions provide the energy 
dissipation which causes damping of stellar pulsations. Ac- 
cordingly, they contribute to the bulk viscosity of matter, 
C. Using the standard definition of the bulk viscosity, the 
energy dissipation rate per unit volume averaged over the 
pulsation period V = 27: /lo can be written as 



i'kin ) = j (div vf ^ ^ 



nbi 

"60 



(5) 



where v is the hydrodynamic velocity associated with 
the pulsations. The latter equality is obtained from con- 
tinuity equation for baryons, nt + riba diyv = (pulsa- 
tions do not change their total number), which yields 
drvv — —rib/ribo = w (jibi/ribo) sinuit. 

The hydrodynamic matter flow implied by the stellar 
pulsations is accompanied by the time variations of the 
local pressure, P{t). The dissipation of the energy of the 
hydrodynamic motion is due to irreversibility of the pe- 
riodic compression-decompression process. Averaged over 
the pulsation period, this dissipation rate in the unit vol- 



^diss 



nb_ 



dtPV. 



(6) 
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For a strictly reversible process, ^£dissy = 0. However, in 
our case the quantities P and V follow variations of Uf, 
in different ways. The specific volume V = l/ni, varies 
instantaneously as ni, varies, i.e., the oscillations of V and 
rib are in phase but the pressure varies with certain phase 
shift. In npeii matter at quasi-equilibrium the pressure can 
be regarded as a function of four variables: ni,, X^, X^, and 
T. Variations of T arc insignificant, for our problem, and 
may be disregarded. Thus, it is sufficient to assume that 
P = P{ni,, Xe, X/^i). Variations of the pressure contain the 
terms oscillating with shifted phases due to the lags of X^ 
and Xf^. 

Let us evaluate the integral We have V — 

— rib/n^Q = (n^i/n^p) sinwi. Thus the only terms in P 
contributing into the energy dissipation are those which 
are proportional to sinojt. At this stage it is convenient 
to use the formalism of complex variables and write 
P = Po+Rc{Pi exp(icjt)}, Xi = Xia+Re{Xii exp(itjt)}, 
where Pq and Xiq are the equilibrium quantities while Pi 
and Xii are small complex amplitudes to be determined. 
We have 



dP 

dx: 



( dP 



\dX, 



X 



^1 1 



(7) 



where all the derivatives are taken at equilibrium. The 
real part of P contains the terms with sin ujt provided the 
amplitudes Xn have imaginary part. 

The change of the lepton fraction Xi is determined 
by the difference of the direct and inverse reaction rates 
given by Eq. (^). The quantity rji in the latter equation 
varies near its equilibrium value ?7;o = as ?/; = + 
Re{77;i exp(ia'i)}, where 



mi 



dill 
drih 



dill 
dXe 



X, 



dm 

dX, 



x^ 



I-lI I 



(8) 



and all the derivatives are again taken at equilibrium. 
Combining the expression n^oXi = F; — f / with Eq. 
and using the formalism of complex variables we obtain 
the two equations Xn = —A; rjn/^iLuni,) (for / = e and /i). 
These two equations supplemented by Eq. (||) constitute 
a system of equations which solution is 



^11 — 



nw {Bii + i ai){Biii' + z a;-) - B^i Bu 



(9) 



where V ^ I and a/ = uuibo/Xi. In analogy with Sawyer 
(|1989D we have introduced the notations: 



Bn 



dXv' 



dm 
dm, 



(10) 



Note that all the derivatives are taken at equilibrium. In 
the absence of muons from Eq. (^) we have = and 
Xei = —{ribi/niio) Cej {Bee + icte)- This is the well known 



limit considered by Sawyer (1989) and Haensel & Schaeffer 

rK . . . 

Generally, Eq. (^ is quite complicated. However, in 

practical applications stellar oscillations are always much 



more frequent than the beta and muon reaction rates 
{\dr]i/dXi\ <C ajnf,o/|A;|) and it is sufficient to use the 
asymptotic form of the solution in the high-frequency limit. 
In this limit the imaginary part of Xn is related to the 
amplitude nti as 



lm{Xii} = -— . 

ribo Lo drib 



(11) 



Combining this equation with that for V (see above) and 
inserting into Eq. (|^) we get the dissipation rate of me- 
chanical energy 



r> \ ^ / '"01 

Cdiss / — \ 



nbi 
nbo 



E 



A; dP dm 

ijj^ dXi dub 



Finally, bearing in mind that (fkin) = ^(^diss) 
Eqs. (H) and (|lj) we have the bulk viscosity 



C = Ce+C, 



dP 



dXi 



dm 

dub' 



(12) 



from 



(13) 



Here we have taken into account that A; and dP/dXi are 
negative and presented the viscosity in the form which 
clearly shows that C; is positive. The expression for 



was obtained by Sawyer (1989). Let us emphasize that 
the viscosity we deal with has meaning of a coefficient in 
the equation that determines the damping rate of stellar 
pulsations averaged over pulsation period (and it cannot 
generally be used in exact hydrodynamical equations of 
fiuid motion). 

Therefore, in the high frequency limit, which is the 
most important in practice, the bulk viscosity C is a sum 
of the partial viscosities Q and produced by the elec- 
tron and muon Urea processes, respectively. This additiv- 
ity rule greatly simplifies evaluation of The values of 
dP/dXi and dm/dub are determined by an equation of 
state as described in Sects. 2.2 and 2.3. The factors A; 
are studied in Sect. 2.4 for the direct Urea reactions. The 
results of analogous consideration for the modified Urea 
reactions will be published elsewhere. 

2.2. Partial bulk viscosity 

Let us discuss briefiy how to calculate the partial bulk 
viscosity Q of npe/j matter for a given equation of state. 
All the quantities in this section and below are essentially 
(quasi)equilibrium values. Thus we will omit the index 0, 
for brevity. 

Since the electrons and muons constitute almost ideal 
gases, the matter energy per baryon can be generally writ- 
ten as 



E = EN{nb,Xp) + XeEe{ne)+Xf,Ef,{nf,), 



(14) 



where Ej^{nb, Xp) is the nucleon energy per baryon, Xp = 
Up /rib is the proton fraction, and Ei{ni) is a lepton energy 
per one lepton (e or /i) . The latter energy is determined by 
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the lepton number density, n; . Owing to charge neutrahty, 
we have Xp = Xp + X^. 

The neutron and proton chemical potentials are given 
by /in = d{rn,Eiq) / dun and /ip — d{ni,Epf)/dnp. The 
derivatives should be taken using Up — Xp rib and n„ — 
(1 — Xp) ni,. This gives fin ^ fJ-p = —dEiq/dXp. The chem- 
ical potentials of electrons or muons arc /i; = {mfc^ + 

p^iC^Y^'^, where ppi — h (Stt^ n;)^^"^ is the Fermi momen- 
tum. Therefore, the difference of chemical potentials is 

dEN{nb,Xp) 



dx. 



Mi- 



lls) 



Now let us calculate C/ from Eq. (10). The derivative 
of fJLi with respect to is evaluated using rt; = X/rif,. The 
result is 



Ci = -Ub 



d'^EM{nb,Xp) 



9 '? 

dub dXp 3 /i; 



(16) 



Using Eq. (|lj) and the standard thermodynamic relations 
we obtain the pressure P = Pn + Pe -\- Pp., where Pm = 
n1 dE^/dnb is the nucleon pressure, while Pe and are 
the well known partial pressures of free gases of e and /i, 
respectively. Direct calculations yields dP/dXi — —UbCi. 
Inserting this derivative into Eq. (|l^) and using the defi- 
nition of C'l we come to a very simple equation 







-7i '^i ■ 



(17) 



Thus, a partial bulk viscosity is expressed through the 
two factors, C; and A;. Calculation of C/ is discussed in 
Sect. 2.3, while A/ is analyzed in Sect. 2.4. 

2.3. Illustrative model of npefi matter 

For illustration, we use a phenomenological equation of 
state proposed by Prakash et al. (1988). According to 



these authors the nucleon energy is presented in the fa- 
miliar form (neglecting small neutron-proton mass differ- 
ence) 



ENirib^Xp) = ENoirib) + S{nb) (1 - 2Xpf 



(18) 



where ENo{nb) — E]y{nb,Xp — 1/2) is the energy of the 
symmetric nuclear matter and S{nb) is the symmetry en- 
ergy. From Eq. (|l^) at equilibrium {rj = 0) we immediately 
obtain /i; = 4(1 — 2Xp)S{nb), and from Eq. (|6|) we have 



Cz = 4 (1 - 2Xp) Tib 



AS 
drift 



(1 - 2Xp) rib -V- 



9 2 

c Pfi 
IJ-i 



drib V 1 — 2Xp J 3 /ii 



2 2 

c Pfi 



(19) 



This is the practical expression for evaluating Ci . The fac- 
tor Ci is not affected by a possible nucleon superfluidity 
which has a negligible effect on the equation of state. The 
relative effect of the superfluidity on the energy per nu- 
cleon is ^ (A//i7v)^ ^ 10~^-10~^, where A is the super- 
fluid energy gap, and /xjv is the nucleon Fermi energy. 



In accordance with Eqs. ( [l9| ) and ([l7|), the bulk viscos- 
ity is determined by the symmetry energy S{nb). At the 
saturation density uq = 0.16 fm^"^ the symmetry energy 
So = S^hq) is measured rather reliably in laboratory (e.g., 
Moeller et al. 1988) but at higher rib it is still unknown. 



Prakash et al. (198?) presented S{nb) in the form: 



S{nb) = UMeY L^/^ _ F{u)\ + SoF{u), 



(20) 



where u — rib/riQ, Sq = 30 MeV, and F{u) satisfies the 
condition F{1) = 1. They proposed three theoretical mod- 
els (I, II and III) for F{u): 



Fi{u)=u, Fu{u) = 



2u' 



1' 



Fiii{u) 



(21) 



and three models for ENQ{nb) in Eq. (|l8|). We do not dis- 
cuss the latter models here because they are not required 
to calculate the particle fractions and the bulk viscosity 



as a function of rib . Following Sawyer (198£) and Haensel 
& Schaeffer ( |l992|) we will use models I and II, for illus- 
tration. Model I gives lower symmetry energy and accord- 
ingly lower excess of neutrons over protons. In contrast to 
the above authors we will allow for appearance of muons. 

The three models for E^olrib) correspond to three dif- 
ferent values of the compression modulus of symmetric 
nuclear matter at saturation, i4ro=120, 180 and 240 MeV. 
If, for instance, we take models I and II of S{nb) and 
the model of EMoirib) with ~ 180 MeV, we have two 
model equations of state of matter (models I and II) in the 
cores of neutron stars. The effective masses of nucleons, 
renormalized by the medium effects, will be set equal to 
0.7 of their bare masses (the same values will be adopted 
in all numerical examples below). The equations of state 
I and II obtained in this way are moderately stiff. The 
maximum neutron star masses for models I and II are 
Mmax = 1.72 M0 and 1.74 M©, respectively. 

The equilibrium fractions of muons and electrons, X^ 
and Xe , can be obtained as numerical solutions of the set 
of the chemical equilibrium equations at given iib'. 



X^ = ]^-AXl'^-X, 



^2/3 „ ^2/3 _ 



? 9 



% (37r^rif, 



,2/3 



0, 



(22) 



where A — hc{3TT^nby^^ / {8 S). If the muons are absent 
(X^ = 0), the second equation should be disregarded, 
while the the first one determines the equilibrium com- 
position of matter 



(23) 



where D = (^'^/27) + (1/16). At given rib the equilibrium 
fraction of electrons in npe/i-matter is always smaller than 
it would be in npe matter, while the equilibrium fraction 
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Fig. 1. Linear response of chemical potential difference 
to baryon number density perturbation, Ci = Ub drji/dnt, 
versus the baryon number density rib for two model equa- 
tions of state proposed by Prakash et al. ( |1988| ) and dis- 
cussed in the text. Solid curves correspond to I — e, while 
dot-and-dash curves to I = ii. Dotted lines present Ce for 
the simplified models of matter without muons. 



of protons is always higher. The threshold of muon ap- 
pearance is determined by the condition /ig = m^c^. For 
models I and II, the muons appear at the baryon number 
density 0.150 fm~'^ and 0.152 fm~'^, respectively. 

In Fig. 1 we plot the factors Ce(nf,) (solid lines) and 
Cfj.{nb) (dot-and-dash lines) which determine the bulk vis- 
cosity for models I and II. The dotted lines show Ce{nb) 
for the simplified models I and II in which appearance of 
muons is artificially forbidden. These results coincide with 
those obtained by Sawyer (1989) and Haensel & Schaeffer 
(1992). They coincide also with the solid lines at den- 
sities below the thresholds of muon appearance but go 
above the solid lines at higher densities (the presence of 
muons affects fractions of electrons and protons). As for 
the factor Cp(nb), it appears in a jump-like manner at the 
muon threshold, initially exceeds Ce{nb) and then tends 
to Ce{nb) with increasing density. 



2.4. Practical expressions for the bulk viscosity of npefi 
matter 

Now let us calculate the factor A;, which enters the bulk 
viscosity (|l^) and determines the asymmetry (^) of the 
rates of direct and inverse reactions of the direct Urea 
processes in npefj, matter. In the absence of superfiuidity, 
the rate of a direct reaction producing a lepton I is given 



by {n 



r, = 



n 



\ (2^)^ 



d^pi d^p^ 



2e,(27r)4 2£,(2^)3 



/l(l-/2)(l-/e) 



(2^)3 6{Ef - E,) 6{Pf - P.) J2 



(24) 



where pj is a nuclcon momentum (j = 1 or 2), p; and 
Si are, respectively, the lepton momentum and energy, 
and El, are the neutrino momentum and energy, S{Ef — Ei) 
and S{Pf~Pi) are the delta functions, which conserve en- 
ergy E and momentum P of the particles in initial (i) and 
final (/) states, |A/p is the squared matrix element of the 
reaction, and fi is an appropriate Fermi-Dirac function, 
/; = {1 exp [{si — fii)/T]}~^ . Equation (p^ ) includes the 
instantaneous chemical potentials fii {i = n^p, I) and does 
not generally require the chemical equilibrium. 

For further analysis we introduce the dimensionless 
quantities: 



T 



(25) 



where the chemical potential difference rji is determined 
by Eq. (||). Thus, the delta function in Eq. ( p^ ) takes the 
form S{Ef ~ Ei) = S{xn — Xp — — + ^), where 
^ = for chemical equilibrium. 

Multidimensional integrals i n Eq. ( p^ ) are standard 
(see, e.g., Shapiro & Teukolsky [1983| ). Equation (|2|) is 
simplified taking into account that nucleons and leptons 
I (e and fi) are strongly degenerate. The main contribu- 
tion into the integral comes from the narrow vicinities of 
momentum space near the Fermi surfaces of these parti- 
cles. The momenta of nucleons and leptons (e or fi) can be 
set equal to their Fermi momenta in all smooth functions. 
The squared matrix element summed over the spin states 
and averaged over orientations of the neutrino momenta 
is 



E 



(26) 



1.436 X 10""*^ erg cm^ 



IS 



Here, G — Gp cos 6*0, Gp = 
the Fermi weak coupling constant, « 1 is the vec- 
tor normalization constant, gA — 1.23 is the axial vec- 
tor normalization constant, and 9c is the Cabibbo angle 
(sin 6*0 — 0.231). Hence the squared matrix element is 
constant and can be taken out of the integral. Further 
procedure consists in the standard energy-momentum de- 
composition of the integration in Eq. (E4). It yields: 

r(=ro/; / = J dxxl J dxndxpdxi f{xn)f{xp)f{xi) 

X S {Xn + Xp + Xi - x^ + Q . (27) 

We have transformed all the blocking factors (1 — f{x)) 
into the Fermi-Dirac functions f{x) by replacing x — > —x. 
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The prefactor Tq is given by (returning to the standard 
physical units) 



47r5 c3 



m'^m*ml{kBT) <dnpi 

1/3 



1.667 X 10 



Here, Tg is temperature in units of 10^ K; to* and 
are, respectively, the effective masses of neutrons and pro- 
tons in dense matter (which differ from the bare nucleon 
masses due to the in-medium effects). Moreover, we have 
defined to* = /ie/c^ ~ PFe/c and to* = [1^,1 <? ~ PFe/c, 
for leptons. The step function equals 1 if the direct 
Urea process is switched on and equals otherwise (Sect. 
1). The direct Urea process of study is switched on if the 
Fermi momenta of the reacting particles satisfy the in- 
equality pFn < (PFp +Pf/) (Lattimer et al. |1991| ). 

It easy to show that the rate F/ of the inverse reaction 
of the direct Urea process (lepton capture) differs from 
the rate of the direct reaction, given by Eq. (p4[), only 
by the argument of the delta function in the expression 
for J (one should replace ^ there). Therefore, the 

difference of the lepton production and capture rates (^ 
can be written as 



Ti -Ti^ To Al, 

nOO 

Jo 

where 

J{x) = JdXndXpdXe f (Xn) f {Xp) f {Xe) 



(29) 
(30) 



(31) 



In a non-superfluid matter, which we consider in this sec- 
tion, the function J{x) is calculated analytically 



J{x) 



2(1 



(32) 



We see that the difference (g^) of the non-equilibrium 
rates of the direct Urea reactions in normal matter is de- 
termined solely by the parameter ^ = rj/T. Moreover, the 
integral (^) is taken analytically for any ^: 



A/ = ll!^e^(0 

60 ^ 



m = 1 



10 



■e + 



1 



■^4. (33) 



This relation, with account for Eqs 
factor A: 



177r2 ^ 177r4 

and (BOh , gives the 



(34) 



In this paper, we do not consider large deviations from 
the chemical equilibrium. We restrict ourselves to the de- 
viations I?]! <C T for which JF(^) w 1. 

Finally, combining Eqs. ( p^ and (jsj), we obtain the 
partial bulk viscosity of npefi matter, Q = (subscript 



26.0 



32 



25.5 



(28) ^ 



p ^ 



g 25.0 
24.5 
24.0 



bjO 



I I I I I I I I I 



I I I I I I I I I 




T=10^K 



no muons — 



4 5 
nb/no 



8 



Fig. 2. The bulk viscosity C^o for models I and II of non- 
superfluid npefi matter (solid lines) induced by the non- 
equilibrium direct Urea processes involving electrons and 
muons versus the baryon number density nt, for T — 10® K 
and u! — 10^ s~^. Dotted lines show the same bulk viscosity 
but for the simplified models of matter without muons. 

refers to non-superfluid matter), induced by a non- 
equilibrium direct Urea process for \t]\ <C T: 



10 



2407r?i^°c3 w2 



^TOpTO, 



ml m: 



^ 8.553 X lO^-* 1121 _ii ( ^ 

no 



1/3 



^ "^9 -"2 



Ci 



100 MeV 



Qnpi gem ^s 



(35) 



where lu^ — u;/(10 s~ ). Figure 2 shows the total bulk 
viscosity = CeO + C/xO of non-superfluid matter versus 
nucleon density Uh. We have used models I and II of npe/i 
matter described in Sect. 2.3. The dotted lines show the 
bulk viscosity for the simplified models in which the muons 
are absent (cf. with Fig. 1). The latter results coincide with 
those reported by Haensel & Schaeffer (1992). 

The results for models I and II are similar. The bulk 
viscosity due to the direct Urea processes is switched on in 
a jump-like manner at the threshold density at which the 
electron direct Urea process becomes operative (this hap- 
pens at nb= 0.414 fm"^ and 0.302 fm"^^ respectively). The 
presence of muons lowers the threshold density (mainly 
due to increasing the number density and Fermi momenta 
of protons). On the other hand, the muons lower the bulk 
viscosity produced by the electron direct Urea process (by 
decreasing Ue)- At larger nt, the total bulk viscosity suf- 
fers the second jump (at rib = 0.503 fm"^ and 0.358 fni^^ 
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for models I and II, respectively). This time it is asso- 
ciated with switching on the muon direct Urea process, 
where muons participate by themselves. The contribution 
of the muon direct Urea into the bulk viscosity is even 
larger than the contribution of the electron direct Urea. 
The total bulk viscosity exceeds the bulk viscosity in npe 
matter. This is natural: the muons introduce additional 
non-equilibrium Urea process which makes stellar matter 
more viscous. 

Finally, let us discuss practical calculation of the bulk 
viscosity. A user possesses usually number densities of dif- 
ferent particles for any given equation of state of npe/i 
matter. This information is sufficient to calculate factors 
Ci numerically from the second equality in Eq. ([l9| ) . Other 
density dependent quantities which enter the expressions 
for the bulk viscosity are also expressed through the num- 
ber densities. Thus, the evaluation of the bulk viscosity 
for any equation of state is not a problem. 

Four our illustrative equations of state the functions 
qi{nb)^{ni/no)^/^ {Ci{nb) /lOOMeY)^, which enter Eq. 
(psj), can be fitted by simple expressions 



Table 1. Studied type of superfluidity 



qi ~ ao + aiu + a2u'^ + 03 u^, 



(36) 



where u = rib/nQ < 15, uq — 0.16 fm ^, and the maximum 
error ^ 0.5%. The fit parameters are: oq — —0.3849, ai = 
0.3338, 02 = 0.0375, as = -0.00098 for / = e and model 
I; ao = -0.3396, oi = 0.4334, 02 = 0.0293, 03 ^ -0.0007 
for I = ji and model I; ao = —1.0491, ai = 1.1788, 02 — 
-0.013, 03 = 0.00044 for Z = e and model II; ao = -0.844, 
ai = 1.2715, a2 = -0.02023, a3 = 0.000675 for / = ^ and 
model II. 

On the other hand, the baryon number number density 
as a function of mass density (for the model equations of 
state with the compression modulus Kq = 180 MeV) can 
be fitted as 



rib = ?^o bip/{po + b2p), 



(37) 



where hi = 1.057 and 62 = 0.0322 for model I; hi = 1.059 
and &2 = 0.0343 for model II. The maximum fit error is 
^ 0.9%, and po = 2.8 x lO^'' g cni^^. These equations 
allow one to calculate the bulk viscosity as a function of 
mass density as required in practical applications. 



3. Bulk viscosity of superfluid matter 

3.1. Superfluid gaps 

Superfluidity of nuclcons in a neutron star core may 
strongly affect the bulk viscosity. Neutrons are believed 
to form Cooper pairs due to their interaction in the 
triplet state, while protons suffer singlet-state pairing 
(Sect. 1). While studying the triplet-state neutron pairing 
one should distinguish the cases of different projections 
mj of nn-pair moment onto a qua ntization axis z (see, 
e.g., Amundsen and Ostgaard 1985 ): \mj\ = 0, 1, 2. The 
actual (energetically most favorable) state of nn-pairs is 



Type 


Pairing state 


FW 


fcBTc/A(0) 


A 
B 
C 


3P2 (mj = 0) 
3P2 (\mj\ = 2) 


1 

(l + Scos^ 1?) 
sin^i? 


0.5669 
0.8416 
0.4926 



not known being extremely sensitive to the (still unknown) 
details of nn interaction. One cannot exclude that this 
state varies with density and is a superposition of states 
with different mj. We will consider the ^P2-state neu- 
tron superfluidity either with mj = or with \mj\ = 2. 
In these two cases the effect of superfluidity on the bulk 
viscosity is qualitatively different. Consideration of the su- 
perfluidity based on mixed nij states is much more com- 
plicated and goes beyond the scope of the present paper. 

Thus we will study three different superfluidity types: 
^So, 3P2 (mj = 0) and i\mj\ = 2) denoted as A, B 
and C, respectively (Table The superfluidity of type A 
may be attributed to any protons, while superfluidity of 
types B and C may be attributed to neutrons. 

Microscopically, superfluidity introduces an energy gap 
d in momentum dependence of the nuclcon energy, e{p). 
Near the Fermi level (|p — Ppl <S^Pf)j this dependence can 
be written as 



fi - + v^{p - pf)'^ at p <Pf, 



e ^ H + yJS'^ +v1{p-pfY at p > Pf 



(38) 



where pp and v-p are the Fermi momentum and Fermi 
velocity of the nucleon, respectively, and /i is the nu- 
cleon chemical potential. In the cases of study one has 

= A2(T)F(i?), where A(T) is the amplitude which de- 
scribes temperature dependence of the gap; F{'d) specifies 
dependence of the gap on the angle d between the particle 
momentum and the z axis (Table ^). In case A the gap is 
isotropic, and S — A(T). In cases B and C the gap depends 
on i). Note that in case C the gap vanishes at the poles of 
the Fermi sphere at any temperature: -Fc(O) = Fc{tt) — 0. 

The gap amplitude A(T) is derived from the stan- 
da rd equ ation of the BCS theory (see, e.g., Yakovlev et 
al. 1999| ). The value of A(0) determines the critical tem- 
perature Tc- The values of fcB?c/A(0) for cases A, B and 
C are given in Table |l|. 

For further analysis it is convenient to introduce the 
dimcnsionless quantities: 



A(T) T S 

V = , T = , V = — . 

r ' Tr T 



(39) 



The dimcnsionless gap y can be presented in the form: 
Ua^va, 2/b = I'd V^l + 3 cos^ yc=«csint?. (40) 
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The dimensionless gap amplitude v depends only on r. In 
case A the quantity v coincides with the isotropic dimen- 
sionless gap, while in cases B and C it represents, respec- 
tively, the minimum and maximum gap (as a function of 
■i?) on the nucleon Fermi surface. The dependence of v on 
T can be fitted as (Levenfish and Yakovlev 1994): 

0.157 1.764' 



T 1.456 



VB 



0.7893 



1.188 



vc 



(2.030 - 0.4903T'* + 0.1727r** 



(41) 



The mean errors of these fits are 1% for all r < 1. 
3. 2. Superfluid reduction factors 

Now let us consider the effects of nucleon superfluidity on 
the bulk viscosity. The dynamics of superfluid is generally 
much more complicated than the dynamics of ordinary 
fluids. Even the motion of matter which consists of par- 
ticles of one species is described by the equations of two- 
fluid hydrodynamics (normal and superfluid components), 
and viscous dissipation of the normal component is deter- 
mined by three coefficients of the second (bulk) viscosity 
(Landau & Lifshitz 1987). Our main assumption is that 
stellar pulsations represent ffuid motion of the first-sound 
type (particularly, temperature variations are neglected) 
in which all constituents of matter move with the same 
hydrodynamical velocity. In this case the hydrodynamical 
equations reduce to the equation of one-fluid hydrodynam- 
ics with one coefficient of the bulk viscosity (C = (2 in the 
notation of Landau & Lifshitz, 1987). 

We will see that superffuidity reduces the bulk viscos- 
ity due to the appearance of energy gaps in the nucleon 
dispersion relation, Eq. (|3^). Quite generally, the bulk vis- 
cosity can be presented in the form 



c 



(42) 



where Qq is a partial bulk viscosity of non-superfluid mat- 
ter, Eq. (^), and Ri is a factor which describes reduction 
of the partial bulk viscosity by superfluidity of nucleons 
1 and 2 involving into a corresponding direct Urea pro- 
cess. If both nucleons, 1 and 2, belong to non-superfluid 
component of matter, we have Ri = 1 and reproduce the 
results of Sect. 2. 

Thus the problem consists in calculating the reduc- 
tion factors Ri . Each factor depends generally on two pa- 
rameters, vi and V2, which are dimensionless gap ampli- 
tudes of nucleons 1 and 2 (and on the type of superffuidity 
of these nucleons). Let us study the effect of superfluid- 
ity on the partial bulk viscosity. For this purpose let us 
reconsider derivation of the bulk viscosity (Sect. 2.1). If 
all constituents of matter have the same macroscopic ve- 
locity, the superffuidity affects noticeably only the factor 



A; in the expression for the bulk viscosity, Eq. (p5|). As 
seen from Eq. (34), the main factor affected by the su- 

which de- 



perffuidity in A; is the integral A/, Eq. (pOD 
scribes the asymmetry of the lepton production and cap- 
ture rates in the direct and inverse reactions of the direct 
Urea process. At ^ <C 1 the integrand of this equation is 
J{xi, - C) ~ Ji^f + C) ~ dJ{x^)/dxi,, where J{xi,) 
is given by Eq. (|3^). Thus, at small deviations from the 
equilibrium one can transform Eq. (BO) to: 



A/o = 4e 



+ 00 



+00 



Axi f{xi) 

dXe f{Xe) S {Xi 



+ CXD 



X2 



dX2 J{X2) 

X,) . (43) 



Here, the index "0" refers to the non-superffuid case, in 
which we have obtained A/q = 177r''^/60. 

Generalization of A; to the superffuid case can be 
achieved by introducing the neutron and proton energy 
gaps into Eq. (^). For convenience, let us deffne the di- 
mensionless quantities 



vp{p-Py) 
T 



£ — /i 

T 



= sign(a;) ^/xPT^ , (44) 
In the absence of superffuid- 



where y is given by Eq. (|3£ 
ity, we have ?/ = and z = x. 

Let the index i = 1 correspond to a nucleon which 
can suffer superfluidity of type A while i — 2 correspond 
to a nucleon which can suffer any superffuidity. A, B or 
C. In order to account for superfluidity in Eq. (^) it is 
sufficient to replace Xi Zi for i = 1 and 2 [in f{xi) 
and in the delta function] and introduce averaging over 
orientations of P2 (analogous procedure is considered in 
detail by Levenflsh & Yakovlev 1994 for the problem of 



superfluid reduction of the neutrino emissivity). Then the 
factor A; can be written as 



Ai 



Azoi?, R{vi,v2) 

f-7r/2 



47r 



J{yi,y2) 



d§ sm^Jiyi,y2) 







where Xio refers to the non-superfluid case, R ■ 
reduction factor in question, and 



(45) 
Ri is the 



J{yi-,V2) = 7 



+ OC 



+c>c 



da;; f{xi) / dxi f{zi) 



+00 



AX2 f{z2) 5{zi + Z2 + xi - x^), (46) 



with 7 = 24G/(177r^). Here, dfi is the solid angle element 
in the direction of pj- 

Thus, we have derived explicit equations 



and (46) 



for calculating the reduction factor R. Calculation is quite 
similar (and in fact, simpler) to that done for the factor 
which describes superfluid reduction of the neutrino emis- 
sivity in the direct Urea process (Levenfish & Yakovlev 
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1994 Yakovlev et al. |999|). The effect of superfluidity on 
the bulk viscosity has also much in common with the ef- 
fect on the emissivity. Thus we omit technical details and 
present only the results and their brief discussion. 

3. 3. Superfluidity of neutrons or protons 

Consider the superfluidity of nucleon of one species, for 
instance, of species 2. In this case R depends on the only 
parameter W2, and we can set zi = xi in Eqs. ( |4^ ) and 
(^). Integration over xi and xi in Eq. (|4^) reduces to 
well-known integrals of the theory of Fermi liquids and 
yields: 



R ^ -f 



tt/2 

dd sin?? 

X [f{z2)B{x^ - 







-t-00 



dXjy Xi/ 



-t-00 



dx-2 



Z2) + f{-Z2)B(x^ + Z2)' 



(47) 



where B{x) = x/(e^ - 1). For r = T2 = T/Tc2 > 1, one 
has R = 1. If superfluidity is strong (r ^ 1, V2 ^ 1), 
the direct Urea process is drastically suppressed by large 
superfluid gap in the nucleon spectrum and reduces the 
bulk viscosity. The asymptotic expressions of i? for r ^ 1 
can be obtained from Eq. (47): 



Ra 



Rc = 



20V2 
177r3-5 

20V3 
5l7r3 
7647r2 
1071 



,3.5 



exp(— u) 



0.221 



exp(- 



^3.5 
0.0367 



exp 



exp 



1.764 

T 

1.188 



4t = 1.7085 

V 



(48) 
(49) 
(50) 



Note that the factors Ra and Rb are suppressed exponen- 
tially with decreasing temperature, whereas Rc varies as 
T^. The latter fact is associated with the presence of gap 



nodes at the Fermi surface (Levenfish and Yakovlev 1994). 

In addition, we have calculated the reduction factors R 
numerically in a wide range of v and propose the expres- 
sions which fit the numerical results (with a mean error of 
^ 1%) and reproduce the asymptotes (|4^)-(|50|): 



Ra 



Rb — 



0.2787 + \/(0.7213)2-f (0.1564 1; 
X exp (2. 9965 - \/(2.9965)2 + i;2 



3.5 



(51) 



0.2854 + v/(0-7146)2-h (0.1418 w)2 
< exp (2.0350 - V(2.0350)2 + ; 



Rc 



0.5-1- (0.1086 1;) 



1 + (0.2347 1')2 + (0.2023 
+ 0.5 exp (1 - ^1 + (0.5i;)2^ 



(53) 



Here, v = va, v = vb and v = wc in the factors Ra, 
Rb and Rc, respectively. Using Eqs. (^) and (|5l|)-(|5^, 
one can easily calculate the reduction factors R for any 
r. These factors are shown in Fig. 3 versus r. We see 



R 



10" 
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10 
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■6 
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Fig. 3. Reduction factors R of the bulk viscosity versus r = 
T/Tc for three superfluidity types A, B and C (Table Q). 

that the reduction can be quite substantial. The strongest 
reduction is provided by superfluidity A and the weakest 
by superfluidity C. For instance, at T = 0.1 Tj, we obtain 
i?A « 2 X 10-^ i?B « 4 X 10--* and i?c » 2 x lO'^. 

3.^. Superfluidity of neutrons and protons 

Let both nucleons, 1 and 2, be superfluid at once, and let 
the superfluidity of nucleon 1 be of type A. In this case R 
can be calculated from Eqs. ( |45| ) and (|4^). Using the delta 
function, we remove the integration over x^, and obtain 

/+00 /•+00 
dxi/(zi)/ dx2fiz2)H{z^+Z2), (54) 
00 •/ — 00 

z -\- X 



H{z) 



dx 



1 



(55) 



Notice that H{z) w z'^/2 as z ^ 00, and H{z) « e^ as 
z —>■ —00. 

First consider the case in which the superfluidities of 
nucleons 1 and 2 are of type A. Using Eqs. ( ^0|) and (45) 
we get 



(52) RAAivi,V2) ^ J{vi,V2) ^ J{V2,VI), 



(56) 



where vi = yi and V2 = 2/2- It is evident that i?AA(0, 0) — 
1. We have also derived the asymptote of Raa in the limit 
of strong superfluidity. Furthermore, we have calculated 
the factor Raa and derived the fit expression which re- 
produces the numerical results and the asymptotes. Both, 
the asymptotes and fits, are given by the complicated ex- 
pressions presented in the Appendix. In Fig. 4 we show 
the curves Raa = const as a function of ti — T/Td and 
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T/T 




T/T 



cp 



Fig. 4. Isolevels of the reduction factor R of the bulk viscosity 
by nucleon superfluidity of type AA versus ri = T/Tci and 
T2 — T/Tc2- In the domain ri > 1, T2 > 1 nucleons 1 and 2 are 
normal and Raa = 1. In the domain ri < 1, r2 > 1 nucleons 
1 are superfluid and nucleons 2 normal, while in the domain 
Ti > 1, T2 < 1 nucleons 1 are normal and 2 superfluid; in these 
domains R depends on one parameter (ri or r2). In the domain 
Ti < 1, r2 < 1 both nucleons 1 and 2 are superfluid at once. 



T2 — T /Tc2- This visualizes the reduction the bulk viscos- 
ity for any T, Tci and Tc2- 

One can observe (Fig. 4) one important property of 
the reduction factor R. If both superfluidities are strong, 

+ r| <C 1, the factor R is mainly determined by the 
larger of the two gaps (by the strongest superfluidity) : 



Ri2 ^ min R2] 



(57) 



Here, Ri and R2 are the reduction factors for the superflu- 
idity of nucleons of one species. The weaker superfluidity 
(with smaller energy gap) produces some additional re- 
duction of the viscosity which is relatively small; this is 
confirmed by the asymptote i?AA given in the Appendix. 
The same effect takes place for the reduction of the neu- 



trino emissivity (e.g., Yakovlev et al. 1999). 

Equation (|4^) can be used also to evaluate R for the 
case in which the nucleons of species 1 suffer pairing of 
type A, while the nucleons of species 2 suffer pairing of 
types B or C. In the Appendix we present the asymptotes 
of the factors i?AB and Rac in the limit of strong super- 
fluidity of both nucleon species (n ^ 1, T2 <C 1). 

The factors i?AB and i?AC can be calculated easily in 
a wide range of ri and T2. The calculation reduces to the 
one-dimensional integration in Eq. (^5|) of the function 
J{yi^y2) fitted by Eq. ([A^). The results are exhibited in 
Figs. 5 and 6. One can see that the dependence of the fac- 
tors i?AB and i?AC on ti and T2 has much in common with 



T/T 




Fig. 5. Same as in Fig. 4 but for the case in which the su- 
perfluidity of nucleons 1 is of type A while the superfluidity of 
nucleons 2 is of type B. 



the dependence of -Raa but i?AB('''i, T2) 7^ ^ab (^2 , ti ) 



and i?Ac('''i: ''2) 7^ -Rac ('''2 , Ti ) . The simple estimate (57) 



turns out to be valid in cases AB and AC as well. However 
since the superfiuidity of type C reduces the factor R in 
a much weaker way than the superfiuidities of types A or 
B, the transition from one dominating superfiuidity to the 
other takes place in a rather wide region of vi and V2 at 
vi ^ \nv2- Accordingly, for V2 ^ vi ^ 1, the reduction 
factor i?Ac exceeds greatly i?AA and i?AB- 

For practical calculations of the bulk viscosity in su- 
perfiuid matter, one needs to know how to evaluate i?AA, 
i?AB and Rac- Corresponding expressions for superfiuid- 
ity of one nucleon species are given in Sect. 3.3. If nucleons 
1 and 2 are superfiuid at once, the reduction factor -Raa 
can be determined easily from the fit equation ([A^). As 
for the reduction factors -Rab and -Rac: we have gener- 
ated their extensive tables. These tables and numerical 
code which generates them are freely distributed. 

Finally, Fig. 7 illustrates reduction of the bulk viscos- 
ity of npefjL matter with decreasing temperature by super- 
fiuidity of neutrons of type B or protons of type A for 
nt, — Auq and oj = 10^ s~^. Thick solid line shows the vis- 
cosity of non-superfiuid matter (cf. with Fig. 2) . Thin solid 
lines exhibit the bulk viscosity suppressed by the proton 
superfluidity at several selected critical temperatures Tcp 
indicated near the curves. The dot-and-dashed line shows 
the effect of neutron superfluidity {Ten = 10^° K) for nor- 
mal protons. We see that the superfluid reduction of the 
bulk viscosity depends drastically on temperature, super- 
fiuidity type, and critical temperatures Ten and Tcp. One 
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T/T 




T/T„ 



cp 



Fig. 6. Same as in Fig. 5 but for the case in which superfluidity 
of nucleons 2 is of type C. 

can hardly expect Ten and Tcp high er tha n 10^" K for rif, as 
large as 4 no (e.g., Yakovlev et al. 1999 ). If so, the super- 
fluid reduction cannot be very large, say, for T ^ 3 x 10^ 
K, but it can reach five orders of magnitude in the case 
of superfluid protons (or six orders of magnitude if n and 
p are superfluid at once, see Fig. 5) for T = 10^ K at 
Ten = Tcp = lO^*' K. It grows exponentially with further 
decrease of T. 




8.0 8.5 9.0 9.5 

Ig T [K] 
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Fig. 7. Bulk viscosity of superfluid npefi matter (model II) 
produced by the electron and muon direct Urea processes at the 
baryon number density — 4 no and lj — as a function 

of temperature for non-superfluid matter (thick solid line) , for 
matter with superfluid protons (solid curves, Tcp = 10^", 10^'^, 
10^ '^ and 10* **) and normal neutrons, and for matter with 
superfluid neutrons (dash-and-dotted curve, Ten ~ 10^*^ K) 
and normal protons. 
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No. 99-02-18099), INTAS (grant No. 96-0542), and KBN (grant 
2 P03D 014 13). 



4. Conclusions 



Appendix 



We have derived practical expressions for the bulk vis- 
cosity of matter in the cores of neutron stars under con- 
ditions in which the bulk viscosity is determined by the 
direct Urea processes. We have paid special attention to 
the case in which dense matter consists of neutrons, pro- 
tons, electrons and muons (Sects. 2.1-2.4). In addition, we 
have studied the reduction of the bulk viscosity by super- 
fluidity of neutrons and protons (Sect. 3). We have an- 
alyzed the cases of singlet-state superfluidity of protons, 
and triplet-state superfluidity of neutrons (without and 
with the nodes of superfluid gaps on the nucleon Fermi 
surface). These cases are most interesting for applications 
(Sect. 1). We have obtained the practical expressions for 
the bulk viscosity of superfluid npefi matter. The results 
can be used for studying the damping of neutron star pul- 
sations (Sect. 1) and in the studies of the gravitational 
radiation driven instabilities in rotating neutron stars. We 
will analyze the bulk viscosity induced by the weaker mod- 
ified Urea processes (||) in a separate paper. 

Acknowledgements. We are grateful to the referee for very use- 
ful critical remarks. Two of the authors (KPL and DGY) ac- 
knowledge hospitality of N. Copernicus Astronomical Center in 



Consider the case in which both nucleon species are super- 
fluid at once. The asymptotic behaviour of the reduction 
factor R given by Eq. (^4|) in the limit of strong superflu- 
idity depends on the gap amplitudes vi and V2 and on the 
type of superfluidities. 

Case AA 

Let both nucleon superfluidities be of type A. If the su- 
perfluidities are strong (wi,W2 ^ 1) and (wi — U2) 3> -^wT 
we obtain: 



Raa 
K 



'..2 , 1 , 2 / «1 



- {vf + 2vi) - -vxvi In 

I \ V2 



(Al) 
(A2) 



where s — \/\v\ — wll- the limiting case vi <C v\ Eq. 
( [A^ ) yields K = vf/60 and reproduces the asymptote 
given by Eq. (|4^). In the opposite hmit y/vi <^ (vi — ^2) ^ 



vi one has K — s /(15w^). I n th e intermediate region 



^2 — ''^il ^ y/v2 the asymptote ( A2) becomes invalid. One 
can show that K ^ y/v^ for vi = W2- 
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We have calculated Raa in a wide range of argu- 
ments vi and V2 and proposed the fit expression valid for 
^/vfTv|<50■. 

u 

= ^.or^o S + D- (A3) 
u + 13.528 

Here, 



^ = 7 (^0 
X exp I 



2.3190 i^i 



0.016689/^2) 



(i) 



1/2 



,1/4 



Vp-q t .0 \ ^« W Vp + Vp-1 



Ki = —Vp~~q, 

6 



P 



~6 



2p = u + 10.397 + Vw^ 



5.4188 u + 7.2987, 



10.397- 



5.4188 u + 7.2987, 



2g = 

2p^ = u + + 3748.6 u + 39.797, 



2pe = M + 0.39417 + \JiiP- + 7.8849 u+ 1.8132, 
D = 0.0080487 {qx qif^ exp(-gi - + 9.9798), 
gi = 2.79205 + ^vl + (2.19785)^, 

52 = 2.79205 



i+ + (2.19785)2 



(A4) 



with u = v\'\'v\ and w = ' 



i?AA given by Eq. (Al), after replacing formally v-i^ yi. 
We have: 



^AB = exp(--ui) 



2 v\v\ 



t 



6^/3 



vj — 6viV2 — 3u; 




{3vl + 11«|) 



2 , / ^^1 +t 



24\/3 



arcsm 



73 



V2 



V2 



(A6) 



4«2. 



where t — \/vf 

The third domain corresponds to ?;i/2 < V2 < vi. 
While averaging over i9 we can split the integral into two 
terms. The first term represents integration over the region 
cos 1^0 < I cosi?| < 1, in which yi < y2, while the second 
term contains integration over the region < |cosi?| < 
cosi?o, in which yi > y2, with cosi?q — {vf — '^2)/(3^'2)- 
The latter term can be taken as an asymptote of the factor 
Rab, since the main contribution into the integral comes, 
again, from | cos7?| <C 1. We have 



-Rab = 7a/ — exp(-wi) 



V2fivi+SV2) 



48^/3 



V2 



(A7) 



Note that the asymptotes given by Eqs. ( |A5|) a nd (|A7[ ) are 



invalid at |w2 — wi| ^ ^/vi- If '^2 ^ ^i, Eq. (A5) reproduces 



Eq. ( '[481) for Ra- For vi <C V2, Eq. ( |A5| ) transforms into 
Eq. (|Aq). Equations ( |A6|) and ( |A7| ) coincide at vi = 2v2- 
If vi — V2, the reduction factor can be estimated as an 



Case AB 

Let nucleons of species 1 suffer pairing of type A and nu- 
cleons of species 2 suffer pairing of type B. In this case 
Vb — y2 = W2\/l + 3cosz?2, see Eq. (^0|). There are three 
domains in the (wi, W2)-plane, in which the asymptotes of 
the reduction factor Rab{vi,V2) in the limit of strong su- 
perfluidity are different. 

The first domain corresponds to V2 > wi, i.e., to y2 > 
yi for all If the both superfluidities are strong {vi ^ 1 
and V2 ^ I) and {v2 — vi) :§> y/v2, the asymptote of Rab 
can be obtained by averaging Eq. ( Al) over 1} after making 
a formal replacement vi — > y2 and V2 ^ vi. Since U2 3> 1 
the main contribution into the integral comes from the 
region in which | cosi9| <g; 1. This allows us to put 7/2 ~ V2 
in all smooth functions under the integral. In this way we 
obtain 



-Rab = 7 



2^/3 



exp(-i;2) 



^ ' ,2 , o„,2x 1 . „,2 , / ^'1 



g (vi + 2«2) - 2^'i^^2 In 



V2 



The second domain corresponds to vi/2 > V2- We have 
Vi > 2/2 for all d in this domain. Then the asymptote of 
Rab is derived by direct averaging over "i? of the asymptote 



Case AC 

Let superfluidity of nucleons 1 be of type A, as before, 
while superfluidity of nucleons 2 be of type C. In this 
case yc — V2 = W2sini?. For deriving Rac in the limit 
of strong superfluidity we will use the asymptote of the 
function H defined by Eq. ( p5|). I f wi > z;2 ^ 1, the main 



contribution into the integral (54) comes from the region, 
in which zi > 22- In this case the asymptote is 



Rac 



7 



■KVl 



exp(-i;i) 



I 2 
y(-i- 



In 



V2 



V7 



6w^wf — J f vi + V2 



2Av2 



(A8) 



If W2 <C wi , this asymptote reproduces the asymptote (|48| ) 
for Ra- For vi — V2, we have Rac ~ jy^Trv/2e~^{3 — 
41n2)w'^/12. Finally, at W2 > I'l ^ 1 the asymptote of 
Rac is 



(A5) i?AC=7J^ exp(-«i)-2 
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